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SUMMARY 


The laminar compressible boundary layer in chemical equilibrium is 
analysed to show the effect of emission of radiation by the boundary layer 
on solutions of the energy equation and on the resulting heat transfer. 
Solutions are obtained at one flight condition, but for several nose 
radii, in a regime where absorption is negligible. The consequent effects 
on heat transfer and bo undary- layer thickness are determined. The concept 
of a separate boundary layer and shock layer is discussed in the light 
of the results obtained. 


INTRODUCTION 


The advent of atmospheric flight at speeds greater than circular 
satellite speed has introduced a number of problems associated with the 
very high temperatures that exist in the air behind the bow shock wave 
of a bluff body. One of these problems concerns the effects of the 
emission of thermal radiation from the chemically dissociated air at 
these high temperatures. From the results of shock tube experiments and 
theory, KIvel and Bailey (ref. l) have made predictions of the emission 
of radiant energy from high-temperature air in equilibrium. Kivel (ref. 
2 ) has applied these results to study stagnation heating by radiation 
emitted from the inviscid isothermal shock layer. 

The Kivel study is a first look at the problem and, of course, 
ignores the effects of gaseous radiation on the structure of the shock 
layer. In Kivel* s work and in other studies of shock layer gaseous 
radiation, the radiation from the boundary layer is generally assumed 
implicitly to be negligible and there is assumed to be no boundary- layer 
interaction with the shock layer. 

Smith (ref. 3 ) in an early analysis, performed from the boundary- 
layer point of view, estimated radiation effects in the boundary layer 
by studying the energy equation in which the convective terms were 
neglected. The results of a numerical example in which water vapor 
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emissivities were used and temperature profries were assumed to be 
unaffected by radiation indicated that radiation emitted by the boundary 
layer has negligible effects on boundary- layer behavior. 

More recently , the problem has been fo: ululated quite elegantly by 
Goulard (ref. 4 ) and again by Tellep and Edwards (ref. 5).. They treat 
the radiation boundary- layer problem in the exact sense by combining the 
astrophysic radiative transfer theory with boundary- layer theory. There 
is some difficulty inherent in solving the resulting set of equations, 
one of which is an integrodifferential equal ion. No solutions were 
attempted in reference 5 - To solve the equations in the stagnation 
region, it is necessary to evaluate the integral part of the integro- 
differential equation over the shock layer. Thus a detailed knowledge 
of the radiating shock layer is essential to solving the boundary layer. 
For this reason, the practicality of trying to solve these equations in 
the stagnation region boundary layer is questionable. 

The purpose of the present paper is to study the radiation problem 
from the boundary- layer point of view, using some approximations that 
make the problem simpler than outlined by reference 5, and yet to use 
recent information concerning the emission properties of air. The aim 
will t>e to verify that there is negligible boundary- layer radiation, that 
temperature profiles are unaltered, that corduction at the wall is 
unaltered, and that wall effects are not felt far from the wall - in 
short, that the shock- layer approach to the radiation problem which 
ignores a boundary- layer effect is satisfactory. In analysing the radiat- 
ing boundary- layer problem, the usual boundary- layer assumptions will 
be made. The results should either provide the verification sought or 
lead to a result that contradicts some of tbe boundary— layer concepts. 
Hopefully, in the event that the contradiction occurs, the results would 
at least suggest a further approach to the problem. 

In this analysis the simplification concerning the radiation effect 
arises if the gas emits but does not absorb radiation to any extent. 

Then the equations to be solved are differential equations instead of 
integrodifferential equations. Thus we set out to solve the boundary- 
layer problem in which there is emission of radiation but no reabsorption. 
Fortunately, this is a very reasonable assumotion for some interesting 
flight conditions. An unpublished study by .C. K. Yoshikawa of Ames 
Research Center concerning the emission and ibsorption of radiation in 
the equilibrium shock layer shows that in a broad flight regime, the 
radiant energy absorbed is negligible compared with that emitted. This 
situation is attributable to the fact that tbe amount of energy absorbed 
depends, in part, on the radiation path leng;h, and for a broad flight 
regime, the path lengths are too short for appreciable absorption. The 
no-absorption regime corresponds roughly to flight at speeds less than 
escape speed at altitudes above 150,000 feet for bodies having a shock 
standoff distance of less than 1 foot. In the present analysis examples 
will be .studied in the flight regime where absorption in the shock layer 

is negligible and it will be assumed that absorption in the boundarv-laver 
air is also negligible. J 
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In addition to the above assumption, it is assumed that the boundary 
layer is in chemical equilibrium, that Prandtl number is constant, and 
that Lewis number is unity. Finally, it is assumed that the only non- 
adiabatic effect in the inviscid shock layer is the emission of radiation. 
Although the analysis is formulated for either two-dimensional or axi- 
symmetric flow, the examples computed will be for axisymmetric flow only 
and will pertain to bodies of different nose radii flying at 31,000 feet 
per second at an altitude of 165,000 feet. For this flight condition, 
it is noted that figure 4 of reference 6 shows the assumption of equilib- 
rium flow in the boundary layer is valid. 


SYMBOLS 

a exponent in equation (A5) 

C Chapman-Rubesin function (eq. (l 9 )) 

c coefficient in radiation emission equation (A5) 

Cp^ specific heat of species i 
Cp frozen specific heat Lcpc Pi 

cq mass fraction of species i 

d exponent in radiation emission equation (A5) 

D coefficient of self diffusion 

Et the rate of total radiant energy emission per unit volume of gas 
f dimensionless stream function (eq. (2l)) 

F dimensionless stream function (eq. (^7)) 

F radiation flux 

g ratio of total enthalpy to total enthalpy at the edge of the 

boundary layer, -JL 
Je 

h static enthalpy 

hp° heat of formation of ith species at 0 ° K 
j total enthalpy 

k thermal conductivity 
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function defined by equation ( 30 ) 

function K evaluated at 7500° K 

heat of dissociation 
c pD 

Levis number, — E. — 
k 

molecular veight of species 

exponent in equation (l); zero for tvc -dimensional flow, unity 
for axisymmetric flov 

exponent in equation (A 5 ) 

pressure 

clu 

Prandtl number, — 

9 k 

heat-transfer rate or energy 

radius of cross section of body 

universal gas constant or body nose raiius 

dimensionless coordinate parallel to tody surface,, equation ( 15 ) 
time 

temperature 

reference temperature (7500° K) 
velocity parallel to body surface 
velocity normal to body surface 

mass rate of production of species i per unit volume 
coordinate parallel to body surface 
coordinate normal to body surface 
intercept in equation (Al) 

velocity gradient at stagnation point (eq. (24) ) 
thermal boundary -layer thickness 

dimensionless coordinate normal to body surface (eq. (l 6 )) 
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coefficient of viscosity- 
gas density 

reference density (sea level air density) 
slope in equation (Al) 
stream function 


JL 

^2 


Superscripts 

derivatives with respect to the independent variable concerned 

Subscripts 

outer edge of the boundary layer 
species i 

reference conditions 

radiant 

stagnation 

conditions in undisturbed stream 

atoms 

molecules 


ANALYSIS 

Basic Differential Equations 


The physical model chosen for analysis is shown in sketch (a) . As 
usual, the region between the bow shock and the body is divided into an 
inviscid shock layer and a boundary layer; both of which emit radiation. 
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The boundary- layer equations expressing continuity of mass, the 
momentum theorem, conservation of energy, and concentration of species 
for a chemically reacting gas are, respectively, 


5 (pur ° n) + ^ (pvr ° m) :: 0 


du 

pU 5x + pv 


du dp 

dy dx 




( 1 ) 

( 2 ) 


pu + pv li = i ) + 

ox oy oy \Pr oy/ oy 


h_ If + A. 


it b s 

M "J5/3y(,2 


Sy 




- div F 


(3) 


and 


dci 


- dci / (3cA 

pu -S + ^ <v pD -s)- *1 


W 


Note that since the boundary layer is assumed to be in chemical 
equilibrium, equation (4) will not be used to determine the degree of 
dissociation; it can be computed directly from local properties. In 
addition, the equation of state of the gas is 


P = 



(5) 


O VJ1 



where 
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M = 




(6) 


and. the definition of total enthalpy is 


d = 



Ur 


Ci hi + %- 


(7) 


where 


h i 



dT + hi° 


( 8 ) 


The exponent m in equation (l) is zero for two-dimensional flow and 
•unity for axisymmetric flow. The third term on the right of equation (3) 
is zero because Le is assumed to be unity. The last term in equation 
(3); div F, represents the rate at which energy is emitted by the unit 
gas volume minus the rate at which energy is absorbed by the unit gas 
volume. It was explained previously that we will be concerned with the 
flight regime in which absorption is negligible compared with emission. 
Then the div F in equation (3) is simply the rate at which energy is 
emitted per unit volume of gas . 


The Radiation Term 


As previously mentioned, reference 1 presents predictions of the 
rate of radiant energy emission per unit volume of equilibrium air as a 
function of temperature and density. In appendix A of the present paper, 
the results of reference 1 are correlated empirically (eq. (A5) ) such 
that the total rate of energy emission per unit volume (which is now 
equated to div F) is given by the expression 


div F = E-t = 



(9) 


where c, d, n, and a are constants appropriately chosen for the regime 
of interest . 
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Boundary Conditions 


The boundary conditions for equations (l) through ( 3 ) are at y = 0; 

u = 0, v = 0, J = (10) 

as y -+ 00 , 

u -» u e ( 11 ) 

Another boundary condition on y is obviously required. For 
solutions of the boundary-layer energy equation without radiation, it is 
usually specified that j approach j e asymptotically as y approaches 
which implies that both the first and second derivatives of j with 
respect to y approach zero at the outer edge of the boundary layer. 
This, of course, is compatible with the usual shock-layer energy equation 
(for an inviscid nonconducting gas without radiation) which is 

p 5t = 0 (12) 

Thus j is constant and there are no gradients of j in the shock 
layer at the outer edge of the boundary layer. 

Now that we have radiation emission from the shock layer and the 
boundary layer, it is necessary to provide a new boundary condition 
compatible with the radiating shock and boundaiy layer. Since the shock 
layer is emitting radiation, we must modify the right-hand side of the 
energy equation (12) . If the emission of radiation is assumed to be the 
only nonadiabatic effect in the shock layer, tie modified energy equation 
is 


P 


au 

Dt 


-div F 


(13) 


or if the Eulerian derivative is written for tie steady state, and the 
boundary- layer coordinate system is assumed to be valid in the shock 
layer near the boundary layer, 



( 1 ^) 


Subsequently, equation (l4) will be transformed in the same way the 
boundary- layer equations are transformed and will be used to provide the 
second boundary condition required. 


< l r\o o\ 
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Transformations 


Equations (l) , ( 2 ) , and (3) and their boundary conditions (lO) and 
(ll) are to be transformed from x and y as independent variables to 
s and t) by means of the Levy transformation (ref. 7), a stream function, 
several definitions and assumptions, and some exterior flow relationships 
as follows: 

The Levy transformation is 


s 


>x 


2m 

p e u e p e r 0 dx 


( 15 ) 


n = 


u e r o 


m 


n/'2sC 



(16) 


A stream function is defined so that 



d\|r _ 

S 7 = 


pvr Q - 


m 


( 17 ) 


and the continuity equation (l) is satisfied. The following quantities 
are defined 


g(n) = -f- 


(18) 


- 3 L m c 

PeM>e 


(19) 


(where C , the Chapman-Rube sin function of reference 8 is assumed 
constant) and 


from which 


f (n) 


U 



f(n) 


i 

■v/ 2 sC 


( 20 ) 


( 21 ) 


In the bluff body stagnation region it is assumed that 


X 


( 22 ) 
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At the outer edge of the boundary layer, the external velocity is 
described by 


u e = px 


( 23 ) 


where from modified Newtonian flow concepts in the stagnation region 




(24) 


It is further assumed that 


PePe ~ (PePe) s 


(25) 


Use of equations (22 ) y 
yields 


(23), and (25) in equations (15) and (l6) 


H(m+i) 

Pp e Pe x 

2(m + l) 


( 26 ) 


A 


n = 


/(m + l)p 

7 c P e Pe 



y 

P <i y 


(27) 


Thus, s is proportional to x 2 and x 4 for two-dimensional and axi- 
symmetric flows, respectively, while q is a function of y weighted 
by the density variation. 

Transformation of equations (2) and (3) to the new independent 
variables s and r\ by means of the defined quantities and assumptions 
results i n_th e following set of ordinary differential equations (for 
constant Pr) 


f 1 1 r + ff " 


ds \ 


1 

P 


) 


e” 

== + fg' 
Pr 



(m + l)Pj e p 



(28) 


(29) 


The right-hand side of equation (28) will be neglected by virtue of the 
qualitative physical argument of reference 9 (based on the fact that the 


VO O VJ1 
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surface temperature is much lower than T e ) . The right-hand side of 
equation (29) will be neglected because Ug 2 < < j e in the stagnation 
region in hypersonic flow. Thus if T and p in equation (29) are 
expressed in terms of g(rj) , similarity solutions of both equations (28) 
and (29) can be obtained. This is done by representing the chemically 
reacting air as a binary mixture of air atoms and air molecules in 
equilibrium. Now letting subscript 1 denote atoms and subscript 2 denote 
molecules, we write the law of mass action in terms of the partial 
pressures (ref. 10, p. 28) 


i. 

5 

) 

) 


V X \~ X = K(T) 

Using the species equation of state 


and the fact that 


in equation (30) yields 


Pi = 



pT 


Mz = 2M X 


-iauu = K (t) 

1 - Cj.2 


Also from reference 10 (p. 30) 


d(ln K) _ M 2 l D 
dT RT 2 


(30) 


(31) 


(32) 


(33) 


(3*0 


The heat of dissociation (noting that the heat of formation of the air 
molecule is zero) 


Z D = hi - h 2 ° = hi° (35) 

is assumed constant. Integrating equation (3^) and substituting the 
result in equation (33) with the boundary conditions K = Kq at T = T 0 
yields 


K(T) = K q 


VdA a 

e r \T 0 y 


= .. > c j . 2 p 
1 - Cl 2 


(36) 


Reference 11 shows that when a mixture of air atoms and air molecules is 
assumed, Cp^ and Cp 2 are both almost constant in a temperature range 
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from 2,000° to 9 .>000° K, and. they differ from one another by less than 
15 percent. Thus we will arbitrarily choose a constant common value of 
C p 1 ~ Cp^ ~ O.318 Btu/lb R between them such that any deviation from 
that common value is less than 10 percent. Then from equation (7) using 
equation (l8) and the fact that vt < < 2j e in the stagnation region 


or 



Cp g T + Cl Z D 

S Q 


= g(n) 


jegCn) - C ! Z P 

c p 2 


(37) 

(38) 


Eliminating c x between equations (36) and (38' yields 


T = 



JeSCn) 



(39) 


Wow T is a function of g (or T}) alone (if we consider that p is 
essentially constant in the stagnation region). By use of equations (5), 
(6) , and (38) , p becomes a function of T and hence, T) 


m 2 P^D 

(*D + JeS - Cp 2 T)RT 


(40) 


If equation (-70) is used in equation (29) , and the right-hand sides of 
equations (28) and (29) are omitted, the differential equations become 


f " ' + ff" = 0 


(41) 


A 


JL 

Pr 


+ fg' - 


r -L_i 

^pZ D I 

_(m + l) P j e P 0 J 

.PqRTUd + j e g - Cp 2 T)_ 


nT a - 


= 0 


(42) 


where T is given by equation (39) . The boundary conditions become 
(from eqs . (10) and (ll)) at t) = 0, 


f = f = 0, 


g 6w 


(43) 


VOOU1 
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at t) -* °°, 


f' -*• 1 (44) 

The outer boundary condition on equation (4-2) is developed by use 
of equation (l4) for the shock layer. Transformation of equation (l4) 
in the same manner as equations (l) , ( 2 ), and ( 3 ) yields the energy 
equation for the shock layer adjacent to the boundary layer 


fg' 


cT 5 - 


_(m + l)Pj e p 0 _ 

|_RT Po (Z d + j e g - Cp 2 T)_ 


(^ 5 ) 


Equation (45) is an expression for the first derivative of g. Applying 
equations (42) and (45) at the edge of the boundary layer (at g = l) 
shows that 


Se" =0 (**6) 

Of course j precisely the same result could have been derived if conduction 
and diffusion effects had been required to vanish at the edge of the 
boundary layer. Equation (4l) is the Blasius equation of reference 12 
if rj is related to the Blasius £ and f ( rj) and its derivations are 
related to the Blasius F(|) and its derivatives by 


4 =-j2 g 
f(n) = fU)//2 
f'(n) = f'U)/2 

f"(T)) - F"(g)/2 J2 






(^7) 


Equations (47) cause the boundary conditions on f ( t)) to be compatible 
with those on F(£) in reference 12. Thus the solution of equation (4l) 
can be obtained from reference 12. It is noted that although a similarity 
type solution of equation (42) is possible, generality is lost in that 
it is necessary to specify a nose radius (for p) and flight condition 
for each solution. 


Heat Transfer to the Wall 


From the solutions of the energy equation, we want to determine the 
heat flux at the wall due to conduction and diffusion and, in addition, 
to determine how much radiant energy originating in the boundary layer 
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strikes the wall (radiant heat transfer to the wall from the region 
between the edge of the boundary layer and the shock is not included in 
this discussion) . To evaluate the radiant energy, we simply assume that 
half the radiant energy emitted from each volume element in the boundary 
layer is directed toward the wall and half toward the shock layer. Then 
in the stagnation region the heat flux at the wall due to conduction, 
diffusion, and radiation emitted from the boundary layer is (for a 
nonreflecting wall) 




k I - pDEhi - U 


,md ( \£\ 

vv 


\nT 


cl- l I dy 


m 


or if equation (48) is rewritten and it is assumed that Le is unity 


q = _ | JL M 

W 6p 


,a 


v 


n P 5 / \ nT’ 

ay 


(49) 


Transforming equation (49) as before and making use of equation (27) 
yields 


A 


lv a - 


Pr 


n 


m + 


l)Cp e q e f3 



\P Z D 

L(Z D + J e S - Cp 2 T) Ep 0 T_ 


lnT 8, -3 


dq 


(50) 


The first term is thus the combined conduction and diffusion heat flux 
and the second te3nn is the radiative heat transfer from boundary— layer 
air for a nonreflecting wall (hereafter referred to as q rg ) . 


Numerical Solution 


Equations (4l) and (42) were solved simultaneously subject to 
boundary conditions (43), (44), and (4 6). The solutions were obtained 
using the Adams Moulton (ref. 13, p. 200) predictor corrector numerical 
integration method programmed for the IBM 704 electronic data processing 
machine. The value of g w ' needed to sta:rt ihe numerical solution of 
equation (42) was determined by trial and errcr so that the boundary 
condition (46) was satisfied. The specified value of g^ corresponds 
to the wall temperature cited below. All examples corresponded to flight 


vo nvn 
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of axisymmetric bodies of varying nose radii at a speed of 31,000 feet 
per second at an altitude of 165,000 feet. The wall temperature was 
fixed at 2000° K. Values of Zp and K 0 were obtained by matching 
equations (36) and (37) with results of reference l4 at a reference 
temperature of 7500° K. 

It became increasingly difficult to obtain solutions as the body 
nose radius increased. The boundary layer thickened, accuracy demanded 
smaller integration steps, and matching boundary conditions became more 
difficult. Solutions were obtained, however, for body nose radii up to 
10 feet. 


DISCUSSION OF THE RESULTS 


In the analysis of this problem, several simplifying assumptions 
have been made. Also, the equation by which the radiation emission is 
calculated is an approximation to experimental results. For these 
reasons, the results of the analysis are not considered to be exact, but 
are intended to show qualitative effects. 

We will examine in some detail the influence of the radiation term 
on the solution of the energy equation. Furthermore , the influence of 
the radiation term on both convective and total heat transfer will be 
evaluated. 

Solutions of the energy equation in terms of total enthalpy are 
shown in figure 1. The upper curve is the enthalpy profile for bodies 
of all radii if there is no radiation term in the energy equation. The 
lower curve corresponds to a solution with radiation emission and pertains 
to a body with a nose radius of 10 feet. Curves of radiation emission 
for smaller nose radii lie between those shown. It is important to note 
that the curve for the no-radiation case goes to unity fairly rapidly, 
while the curve with radiation approaches unity very slowly (at values 
of 7} greater than those shown) . The emission of radiation effectively 
thickens the thermal boundary layer (thickness here being defined by 
that point where equation (47) is satisfied) thus including more of the 
shock layer air within the boundary layer. 

The temperature profiles corresponding to the enthalpy profiles are 
shown in figure 2. Again, the upper curve corresponds to the usual 
solution of equation (43) for all body sizes if the radiation term is 
omitted. It exhibits the usual inflections for equilibrium boundary 
layers (these inflections are also evident in similar plots presented In 
ref . 15) . Again, the upper curve for the no-radiation case converges to 
unity very rapidly. The lower curves in figure 2 correspond to solutions 
of equation (43) with the radiation term included and are for nose radii 
of 3 and. 10 feet. Here the difference between the radiation and no- 
radiation solutions is somewhat more pronounced than the differences 
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observable on the previous figure. The thermal boundary layer is thicker 
and, of course , conduction and diffusion effects are felt much farther 
from the vail. 

It should also be noted that the emissior of radiation by the boundary 
layer diminishes the temperature and enthalpy gradients at the vail small 
amounts (actually somevhat too small to be seen on fig. 2) . Such a 
result is to be expected since the greater proportion of radiation is 
emitted from the hotter air not in close proximity to the vail. It can 
be anticipated that the corresponding reducticn in convective heat 
transfer Is more than offset by the added radiant heating load. More 
vill be said of this in connection vith a subsequent figure. 

Figure 3 shovs the distribution of radiant energy emitted in the 
boundary layer normalized to the value at the edge of the boundary layer 
for the nose radii of 5> and 10 feet. The figure shovs that the 
emission rate rises rapidly up to r) ^ 5, after vhich it rises more 
gradually to unity at the boundary- layer edge. If tj vere linearly 
related to y, the area under a given curve up to any value of r\ vould 
be proportional to ( l T ( r i) > the rate of total radiant energy emission 
from the layer of air betveen the vail and T], Hovever, because of the 
compressibility effect, t) is not linearly related to y (see eqs. (l6) 
and (27)) , and to obtain q (tj) as a function of T), it necessary to 
perform the integration shown in equation (51), vhere q r (T])/q r is the 

ratio of the radiant energy emitted, from the layer of air betveen the 
vail and r\ to the total radiant energy emitted by the boundary layer: 


q. r (n) 


T U 


Mgp^p 

l_(*D + J e S - Cp 2 T) P:) RT 


nnT 8 --! 


dr) 


f % Ti 


L 

_(Zp + j e g - Cp 2 T) p 3 RT_ 


inT a -1 


dt] 


(51) 


The ratio is plotted in figure 4. Without the radiation term in the 
energy equation, the solution converged at a value of q of roughly 10 
for all nose radii. Thus the figure shows that most of the radiation 
emitted by the boundary layer is from air that for the no-radiation case 
would have been outside the boundary layer. 

The ratio of the total heating load (computed from eq. (51)) including 
incident radiation from the boundary layer to neat transfer without 
radiation is plotted as a function of nose radius in figure 5. It should 
be noted that the comparison does not include the radiation contribution 
from the region between the edge of the boundary layer and the shock wave. 
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It was mentioned previously that the slight reduction of convective 
heating due to boundary- layer radiation emission is more than offset by 
the added radiative heating. This is shown in the figure by the fact 
that the ratio plotted is always greater than unity. For the 3-foot 
nose radius, the total heat transfer is increased by almost one third 
because of radiation heating, and for the 10-foot nose radius, it is 
more than doubled. 

Now that the detailed solutions of the boundary— layer energy equation 
in the no-absorption regime have been discussed, it is pertinent that 
an evaluation be made of the boundary- layer approach for studying gas- 
dynamic effects in the stagnation region when self-radiation is important. 
In the analysis, several simplifying assumptions and concepts have been 
used. It is not usually possible to assess the full significance of all 
the assumptions in advance, but one hopes to be able to tell by the 
resulting solutions just how reasonable his assumptions were. In the 
present case , one further step has been made in the usual boundary -layer 
analysis in order to include the effect of radiation emission. In the 
resulting solutions the momentum boundary layer was unchanged. However, 
the solutions of the energy equation showed that its boundary conditions 
were satisfied at distances from the wall several times that of the 
corresponding no-radiation case. Indeed, the thermal boundary- layer 
thicknesses were comparable to the thicknesses usually calculated for the 
inviscid adiabatic shock layer. For this reason, the concept of a 
thermal boundary layer (to which conduction and diffusion effects are 
limited) separate and distinct from a shock layer loses its usef uln ess - 

The results leading to this conclusion were obtained for the no- 
absorption case. An attempt was also made to solve the very strong 
absorption case. The results of that effort will not be discussed 
except to say that the remarks of the preceding paragraph also apply to 
the very strong absorption case , and additionally that the convective 
heating at the wall is greatly altered by radiation effects for the very 
strong absorption case. 

The next approach to studying the gaseous radiation problem for 
the bluff body suggests itself. It is, of course, to solve the flow field 
from the wall to the shock including viscous, conduction, diffusion, and 
radiation effects in the flow equations. It is anticipated that such 
solutions would also give the shock standoff distance. Thus an attempt 
to obtain solutions from the body surface to the shock seems to be in 
order. In that light, the present study may be regarded as a first step 
toward this goal. 


CONCLUDING REMARKS 


The effects of radiation emission from the boundary- layer air have 
been studied by including an emission term in the energy equation. This 
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term was expressed by a relationship develoj ed from experimental results 
of Kivel and Bailey. An appropriate boundary condition on the energy 
equation has been derived which limits condi.ction and diffusion effects 
to the boundary layer, and which allows the shock layer total enthalpy 
to vary as a result of the emission of radiation from the shock layer. 

For the flight condition used in the sclutions the emission of 
radiant energy diminishes the total enthalpy by a small amount. Similarly, 
the enthalpy gradient at the wall and thus Ihe convective heat transfer 
are diminished slightly. However, the reduction in convective heat 
transfer is more than offset by an added radiant heat-transfer load. 

Existence of an enthalpy gradient at the outer edge of the boundary 
layer and the self-emission of radiant energy increase the thickness of 
the thermal boundary layer by a large amount. More air (which is at an 
elevated temperature and is emitting radiation) is associated with the 
boundary layer. Conduction and diffusion effects are felt farther from 
the wall. For this reason, the concept of e thermal boundary layer 
separate and distinct from the shock layer vhen self-radiation is consid- 
ered loses its usefulness in gaseous radiation studies that are more 
refined than those now in existence. At present the only way to calculate 
the gaseous radiation contribution to aerodynamic heating is to neglect 
the boundary layer and surface effects and use the shock layer results 
alone. Such an approach to estimating the radiative heat transfer, of 
course, gives no information of radiative effects on flow field structure, 
standoff distance or convective heating. Although the convective heating 
is not significantly affected in the present analysis for the no-absorption 
flight regime, it is greatly altered in the very strong absorption regime, 
and other methods must be used to study the problem. It is anticipated 
that in the intermediate regime (where self-absorption is significant 
but not excessive) , convective heat transfer is substantially altered by 
the radiation effects, and that boundary layer and shock layer are again 
not separable. It is suggested that the next approach be to solve the 
flow field from the body to the shock without imposing a boundary- layer 
concept . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., June 20, 195l 
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appendix a 

CORRELATION OF RADIANT EMISSION DATA 


Predictions of the total rate of radiant energy emission per unit 
volume from high-temperature air are given in reference 1. To utilize 
this information in the present analysis, it was necessary to derive an 
analytic expression for the variation of the radiant energy with temper' 
ature and density. Results from reference 1 are shown on figure 6 as 
dashed lines. Each curve (which was determined by five points) was 
replaced by its individual least squares best fit straight line. These 
are the straight lines shown in figure 6 and are represented by the 
expression 


logic (jf) = lo glo a(T) + T(T)lo glo (jf) (Al) 

or 

If (P,T) = a(T)(^) T(T) ( A2 ) 

where a(T) and t(t) are listed on the figure. Secondly, all the slopes 
and all the intercepts logio T (T) and log 10 a(T) were fitted by the least 
squares best fit straight line; that is, the lines 


l°gio T (T) = logion + a log 10 T (A 3 ) 


and 


logioa(T) 


logio 



+ d log 10 T 


(A4) 


were obtained. Combining equations (A2) , (A 3 ), and (A4) yields 


Et(p,T) 


cT 1 



(A5) 


Thus all the lines in figure 7 were obtained from the single expression 
(A 5 ) where c, d, n, and a are shown in the figure (where T is in 
°K, and Et is watts/cm 3 ) . 
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Although the representation of the prediction of reference 1 hy 
equation (A5) is fair in figure 7> it can he improved in a restricted 
range of temperatures and densities if the technique is applied to the 
data of that range alone. The results of using the data in the temper- 
ature range 1000° to 8000° K and the density range 10" 3 < p/p Q < 10" 1 
are shovn in figure 8. Equation (A5) vith constants shovn on figure 8 
was used in the calculations of this paper. 


VO O VJ1 > 
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Figure 2.- Temperature profiles. 
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Figure 5*- Ratio of total boundary- layer heat transfer to wall with 
radiation to that without radiation as a function of nose radius 
(flight speed 31,000 ft/sec at 165,000 ft altitude, T w = 2000° K] 
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Figure 7*~ General correlation of all data. 











